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Abstract 

We investigate the interaction between P baryons in the ^S'o channel from 2-1-1 flavor lattice QCD 
simulations. On the basis of the HAL QCD method, the PP potential is extracted from the Nambu- 
Bethe-Salpeter wave function calculated on the lattice by using the PACS-CS gauge configurations with 
the lattice spacing a ~ 0.09 fm, the lattice volume L ~ 2.9 fm and the quark masses corresponding 
to niTr ~ 700 MeV and mn — 1970 MeV. The PP potential has a repulsive core at short distance and 
an attractive well at intermediate distance. Accordingly, the phase shift obtained from the potential 
shows moderate attraction at low energies. Our data indicate that the PP system with the present 
quark masses may appear close to the unitary limit where the scattering length diverges. 

Introduction Strange dibaryons have been attracting considerable interests both theoretically and 
experimentally in hadron physics. In particular, the 77-dibaryon with (strangeness)=—2 [T] and the AQ 
dibaryon with (strangeness)=—3 [2] are considered to be the promising dibaryon states due to the absence 
of Pauli repulsions among valence quarks at short distance (see the reviews, mm- In recent years, the 
numerical and theoretical progresses in lattice gauge theories made it possible to attack such a problem 
directly from quantum chromodynamics (QCD) (See e.g.[5l [6l [7] and references therein.) 

The purpose of this letter is to extend our previous works on the (strangeness)=—2 systems such 
as 77-dibaryon [5] and NE [8] and the (strangeness)=—3 system as NQ [9] to the (strangeness)=—6 QQ 
system in 2-|-l flavor lattice QCD. In our approach (the HAL QCD method), the baryon-baryon potential is 
extracted from the Nambu-Bethe-Salpeter (NBS) wave function calculated on the lattice: Such a potential 
deduced in lattice QCD is guaranteed to reproduce physical observables (e.g. the scattering phase shift) 
by construction m- The HAL QCD method has several advantages over the conventional finite volume 
method HI]: (i) The finite volume effect is highly suppressed due to the short range character of baryon 
potentials, (ii) the ground state saturation of the two-particle system is not required for extracting the 
potential, since the same potential distates all the scattering states on the lattice, and (iii) physics behind 
the baryon-baryon interaction can be easily grasped from the spatial structure of the potential. Further 
details on these points are discussed in [Ellll]. 
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There exit several investigations so far on the interaction using the phenomenological quark models: 
Some studies show strong attraction which may cause a QQ bound state HaE!, while other studies show 
weak repulsion [Ml [I7]. A recent lattice QCD analysis of the 00 scattering length m by using the 
standard finite volume method m shows weak repulsion in the jS-wave scattering with the scattering 
length oqq = —0.16 ± 0.22fm: This indicates that the 00 system is unlikely to have a strongly bound 
dibaryon, although the large error prevents us to make a firm conclusion about details of the interaction. 

The HAL QCD potential Let us first recapitulate the essential part of the HAL QCD method to be 
used for extracting the 00 potential. The basic quantity is the equal-time NBS wave function with the 
Euclidean time t, 

^a,k{i^mpAt,o) |oo,iy), ( 1 ) 

where |00,kL) is an exact (strangeness) = —6 state with zero total momentum. The total energy of the 
system is given hy W = 2^J with the 0 baryon mass mo and the relative momentum p. Local 
interpolating operators for the 0 baryon, 0 (x) and 0 (x), are taken to be 

^a,k{x) = Sa,a{x)sb{x)'ykiCs'^{x), ( 2 ) 

where a, b and c are color indices, is the totally anti-symmetric tensor, 7 ^ represents the Dirac matrix, 
and a is the spinor index. The charge conjugation matrix is taken as C = 7472 = —C~^ = —C'^ = —C^ 
in the Euclidean space. 

An important property of the NBS wave function is its asymptotic behavior at large distance, denoted 
simply as 

pr 

where p = \p\, r = |f^, D is the solid angle of r, Ylm is the spherical harmonic function, and L is the 
orbital angular momentum. Although Eq.([3|) looks like a quantum mechanical formula, it can be derived 
from the unitarity of the S'-matrix in quantum field theory with 5l{p) being the scattering phase shift for 
given quantum numbers in QCD [10] . 

Our next task is to define the potential from which this scattering phase shift can be calculated. In 
the HAL QCD method, such a potential is defined through the Schrodinger type equation, 

{E-Ho)A^{r) = J (4) 

where Ho = — 2 ^^^ is the free Hamiltonian, = ra^ij^ is the reduced mass, and E = i® 

kinetic energy. Here U (r, r^) is the non-local but energy independent potential, which can be expanded 
in terms of the non-locality (the velocity or derivative expansion) [19]. The convergence of the velocity 
expansion at low energies has been investigated previously for the nucleon-nucleon scattering m and 
the pion-pion scattering [ 21 ]. Since we consider the low-energy S'-wave scattering much below the meson 
production threshold in this letter, we take only the leading-order local potential, V{r), in the expansion, 
U{r,E) = Vir)5{r-T)+0{V). 

Note that V (r) is an effective central potential, which contains not only the genuine central part but 
also the tensor part implicitly m- 

The NBS wave function can be extracted from the asymptotic temporal behavior of the four-point 
(4-pt) function, 

F(f, t - to) = (0|L!(t, r)n{t, 0) J(to)|0) = J^(0|L!(t, f)L!(t, 0)|L!L!, 1T„)(L!D, Wn\J{to)\0) + ■■■ 
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{t — tQ ^ oo) 


(5) 



+ • • • ~ aoV''^°(r)e 


where i 7 (io) is the wall source operator which creates the flfl state at to, an is the matrix element defined 
by (no, VF„|^7(0)|0), Wn are discrete QCD eigen-energies in the finite volume below inelastic threshold, 
and Wq is the lowest eigen-energy. The ellipses in the above equation represent inelastic contributions in 
the on system, which are suppressed for large t — to. 

Eq. ([5]) shows that the NBS wave function for the ground state can be extracted from the 4-pt function 
at large t, as long as other contributions from Wn>i > Wo can be neglected. In practice, however, the 
increasing statistical errors of the 4-pt function at large t make it difficult to achieve the ground state 
saturation with reasonable accuracy [22]. Moreover, as the volume increases, larger and larger t becomes 
necessary to extract These problems can be simultaneously avoided by the time-dependent HAL 

QCD method introduced in m, where we start with the so-called i 2 -correlator. 






( 6 ) 


^ AIT^ 

with AH^i = Wn — ^niQ. Since AITn = — - -we have 

mn 4mf2 




R{f,t) ~ = 


fn 1 9^ 

mn 4mQ dt"^ 




)an 0 


= ( -- 
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mn Amn dt'^ 


)R{f,t) + / drr U{r,r)R{f,t), 


(7) 


where we have replaced hy U — V'^/tuq using Eq.Q. We then arrive at the time-dependent 

equation, 

m ^ I dh'u(f.f)mf.t), (8) 

where implies that we have neglected inelastic contributions by taking sufficiently large t — to- Eq.([5D 
gives U (r, H) directly from F{f, t). This method no more requires the ground state saturation, so that the 
moderately large values of t — tg which suppress inelastic contributions are enough for a reliable extraction 
of the potential. Then, in the leading order of the velocity expansion, we obtain 


H(f) 


1 

mn 


- 




R{f, t) 


(9) 


Interpolating operators for QQ system The present system can be characterized by the total 
spin (S'), the orbital angular momentum (L), the total angular momentum (J), and the parity P. The 
asymptotic HD state with given L and S has a factor (— 1 )'^+^+^ under the exchange of two H’s, so that 
S -|- L must be even due to the Fermi statistics of H baryons. In Tabled] we show low-J channels 
which appear for given conserved quantum numbers J and P. In this letter, we employ the wall source, 
L = 0 with S = 0 at t = to) which creates the = O'*" state, so that only the upper left corner of this 
table is considered. Then, both ^Sq and ^Dq channels appear after the QCD interactions at t > to- As 
mentioned before, we determine only the effective central potential from the ^Sq channel at t > to) where 
effects of the ^Dq state are implicitly included. 

The interpolating operators for Q with 5 = 3/2 and 5^ = ±3/2, ±1/2 read 

Da 3 = — {tpT 1 , ( 10 ) 

2 ’ 2 2 
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P = + 

P=- 

o 

II 

'So, 

"To, ^To 

J = 1 

"Pi 

"Ti, ^Ti 

J = 2 

"S2, 'P2, ®P2, ®G2 

"T2, ^P2, "T2, ^T2 , ^T2 

J = 3 

"P3, "G3 

"P3, "T3, ^T3 , ^T3 

J = A 

"P4, 'G4, "G4, "/4 

^P4, "T4, "T4, ^T4, ^4 


Table 1: A relation between conserved quantum numbers (J and P) and quantum numbers in the asymp¬ 
totic nn channel. 


Ds 1 = ^[\/2(T/’r^i/>)V'i-h ('i/’r+V’)V’_i]) 

2’2 .^3 2 2 

(11) 

Ds 1 = 1-1-(V’r-V')V’i]) 

2’ 2 V3 2 2 

(12) 

Ds 3 = (V’r_V')i/’ 1, 

2 ’ 2 2 

(13) 


where r± = ± iC'y^) and Tz = We take only upper two components in the Dirac 

representation for the quark operators, so that the D operator does not have the S' = 1/2 component. 
Combining these operators, the interpolating operators for the DD system with the total spin S = 3, 2,1,0 
with Sz = 0 are given by 


3 3 Da _3 -|- 3D 3 1 Da _1 3D 3 _iD3 1 -l-Da _3D3 3 ), 

2’2 2 ’ 2 2’2 2 ’ 2 2 ’ 2 2‘2 2 ’ 2 2’2 


(DD)2 0 — ~(D33D3_3-|-D3]_D3_1 — D3_1.D31 — D3_3D33), 

’ 2 2’2 2 ’ 2 2’2 2 ’ 2 2 ’ 2 2'2 2 ’ 2 2’2 


3 3 D 3 _3 — Ds iDa _1 — Ds _iD3 1 3D 3 _3D3 3 ), 

2’2 2 ’ 2 2’2 2 ’ 2 2 ’ 2 2'2 2 ’ 2 2’2 


(DD)n n = —(D 3 3 D 3 3 —Da iDa i +D 3 1 D 3 1 —Da 3 D 3 3 ). 

’ 2 2’2 2 ' 2 2’2 2 ' 2 2 ’ 2 2'2 2 ’ 2 2'2 



(DD)3,o - 


(14) 

(15) 

(16) 
(17) 


In this letter, we use Eg. lfTTl) to calculate the S = 0, Sz = 0 state. 

To extract the L = 0 state at the sink t on the lattice, we employ the cubic-group projection defined 


by 


da 


9 


9 


i 


(18) 


where a represents an irreducible representation of the cubic group with dimension da, Ri is an element 
of the cubic group acting on r of D operators, D°‘{Ri) is the corresponding matrix in the irreducible 
representation acting on the spin of D operators, and g is the order of the cubic group. We use the Ai 
projection, which generates the L = 0 as well as the L = 4, 6, • • • states, where the L = 4,6, • • • components 
are expected to be negligibly small. For example, we have 


P/'iDo,o(r1 = ^f;Do,o(i?.rl. 

i=l 


(19) 


The DD potential We employ 399 gauge configurations generated by the PACS-CS Collaboration with 
the renormalization group improved gauge action and the non-perturbatively 0{a) improved Wilson quark 
action in 2-|-l flavor QCD [23]. These configurations were obtained at /3 = 1.90 (a = 0.0907(13) fm) on 
the 32^ X 64 lattice, whose physical extension is L = 2.9 fm, with the hopping parameters Kud = 0.13700 
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Figure 1: The i?-correlator in the ^Sq channel as a function of the relative distance r for t — to = 12. 


and Kg = 0.13640, corresponding to m-^ = 701(5) MeV and mn = 1966(6) MeV. We employ the wall 
quark-source with Coulomb gauge fixing. The periodic (Dirichlet) boundary condition is imposed in 
spacial (temporal) direction. To improve the statistics, we perform the measurement at 64 source time 
slices for each configuration, where the unified contraction algorithm [23] is used to calculate the NBS 
wave functions. Statistical errors are estimated by the Jack-Knife method. We make analyses with the 
bin sizes of 1, 3, 7, 19, 21 and 57, and the bin size dependence is found to be negligible. Hereafter, we 
show the results obtained with the bin size of 1, unless otherwise indicated. 

Fig. [1] shows the ii-correlator (Eq.([6|)) in the channel as a function of r for t — to = 12, where 
the i?-correlator is normalized to be 1 at r = 2.5 fm. We find that the i?-correlator is enhanced at 
intermediate distance, while it is suppressed at short distance. The latter behavior is consistent with the 
partial Pauli blocking in the quark level, similar to the situations in the nucleon-nucleon force [El- 

Shown in Fig. [2] is the effective central potential Vc{r) between H baryons at t — to = 12 in the ^So 
channel. The Laplacian term and the time derivative term calculated from the i?-correlator in the right 
hand side of Eq. ([8|) are separately plotted in the figure, together with the total potential. We here 
approximate the time derivative term in Eq.Q by ■§iR{t) = i) = R{t + 1) + R{t — 

1) — 2R{t). We hnd that the time derivative terms have sizable contributions to the total potential: The 
1st derivative in t dominates over the 2nd derivative in t. The latter corresponding to the relativistic 
effect is negligible. Eig. [3]shows the time dependence of Vc{r) at t — to = 11, 12, 13. This particular region 
of t is chosen to suppress contaminations from excited states in the single H propagator at smaller t and 
simultaneously to avoid statistical errors at larger t. We observe that the potential is nearly independent 
of t within statistical errors, as expected in the time-dependent method |12j . 

The effective central potential 14 (r) in Eig. [3] has a repulsive core at short distance and an attractive 
well at intermediate distance. This is qualitatively similar to the NN and HE interactions, where the 
partial Pauli blocking in the quark level takes place and the i?-correlators are suppressed at short distance. 
On the other hand, no Pauli blocking is active for the H dibaryon and the NO, dibaryon, so that there is 
no repulsive core in these cases UM- 


5 





















2000 


1500 


> 

(D 

J 1000 
> 

500 


0 


0 0.5 1 1.5 2 2.5 

r[fm] 

Figure 2: The effective central potential for the flfl system in the ^Sq channel at t —to = 12. We separately 
plot the Laplacian term (blue), the time derivative term (green) and the total (red). 
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The nn phase shift To obtain the 5-wave phase shift, we fit the potential in FigjS] using a function 
which contains two Gaussian terms plus the Yukawa squared term with a form factor m, given by 


o-ajr 


V,(r) = + a 5 (l - 

r 


lim V(r) = ai + a^. 

r —>-0 


( 20 ) 


This 2 Gauss -|- (Yukawa)^ form gives a fairly good fit with y^/d.o.f ~ 0.50 at t — to = 12- The re¬ 
sulting parameters are oi = 1.69(6) x 10^ MeV, 02 = 1.24(3) x 10^fm“^, 03 = 4.44(68) x lO^MeV, 
04 = 5.68(1.31) fm-^ 05 = -7.06(14.64) x lO^MeV, oo = 6.25(5.77) x 10"^ Im'^ 07 = 3.43(30) fm"! at 
t-to = 12 . 

Using the fitted potential, we solve the Schrodinger equation in the infinite volume to calculate the 
scattering phase shift difk) of the HU system in ^5o channel by the formula with L = 0 , 


tan5i(/c) 


V’fc(a^2) sin(/cxi - ^) - s\n{kx2 - %) 

hm - f - f —, 

X1,X2^-00 ll)^[xi) COs{kX2 - - 'Ipk{x2) COs(/cXl - 


( 21 ) 


where tpk is the wave function and k is the magnitude of the momentum. 

Fig. IHshows the phase shift as a function of the kinetic energy, E = k'^/mn. The result indicates that 
the nil interaction is attractive at low energies, while the existence of the bound state is inconclusive 
because of large statistical errors near A; = 0. Indeed, the central value of (5o(^ = 0) at t — to = H is 
180°, while it becomes zero at t — to = 12,13. Due to the same reason, the scattering length and the 
effective range cannot be extracted reliably from this phase shift. A possible physical interpretation of 
such a situation is that the DD system at the present quark masses may appear close to the unitary limit 
where the scattering length diverges and changes its sign. 
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Figure 3: The effective central potential 14(r) in the ^S'o channel at t —to = H (blue), 12 (red), 13 (green). 



Conclusion In this letter, we have calculated the effective central potential and the scattering phase 
shift between Q baryons in the channel by using the HAL QCD method. We have found that the 
fill potential has short range repulsion and intermediate range attraction just like the nucleon-nucleon 
potential. The short range repulsion of this system is a reflection of the partial Pauli blocking in the 
quark level similar to the nucleon-nucleon potential, and is in contrast to the cases of the LT-dibaryon or 
Nil system with no repulsion. The HH interaction is attractive at low energies, but is not strong enough 
to form a tightly bound dibaryon at quark masses corresponding to — 700 MeV and niQ ~ 1970 MeV. 
Rather, the system may appear close to the unitary limit at these quark masses. We plan to carry out 
the 2-|-l flavor simulations at the physical quark masses, in order to investigate whether the attraction 
found in the present study increases or decreases toward the physical quark masses. 
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